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ABSTRACT 

Aims. The aim of this study is to investigate the nature of radio sources associated with young stellar objects (YSOs) belonging to the 
R CrA cluster. By combining the centimetre radio data with the wealth of shorter wavelength data accumulated recently we wish to 
refine estimates of the evolutionary stages of the YSOs. 

Methods. The region was imaged at 3, 6, and 20 cm using the Australia Telescope Compact Array. Fluxes and spectral indices for the 
brightest radio sources were derived from these observations. Some of the 3 and 6 cm data were obtained simultaneously, which is 
useful for reliable spectral index determinations of variable sources. 

Results. Eight pointlike sources were detected. Seven of them can be assigned to YSOs, which have counterparts in the X-rays, 
infrared or submm. One of the YSOs, Radio Source 9, is a Class candidate, and another, IRS 7B, is suggested to be in the Class 
0/T transition stage. IRS 7B is associated with extended radio lobes at 6 and 20 cm. The lobes may have a gyrosynchrotron emission 
component, which could be understood in terms of Fermi accleration in shocks. The Class I objects detected here seem to be a mixed 
lot. One of these, the wide binary IRS 5, shows a negative spectral index, rapid variability, and a high degree of circular polarisation 
with V/I w 33 % on one of the days of observation. These signs of magnetic activity suggest that at least one of the binary components 
has advanced beyond the Class I stage. The radio source without YSO assigment. Radio Source 5, has been suggested to be a brown 
dwarf. The radio properties, in particular its persistent strong emission, do not support this classification. 

Conclusions. The radio characteristics of the detected YSOs roughly agree with the scheme where the dominant emission mechanism 
changes with age. The heterogeneity of the Class I group can possibly be explained by a drastic decline in the jet activity during this 
stage, which also changes the efficiency of free-free absorption around the protostar 
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1. Introduction 

The presence of a protostar in a molecular cloud core is indi- 
cated by compact centimetre radio emission, a bipolar molec- 
ular outflow, and/or near- or mid-infrared (NIR/MIR) emis- 
sion (lAndre et al. 2000l l. Centimetre continuum radiation ob- 
served toward YSOs originates in most cases from thermal 
free-free emission from partially ionized stellar winds or jets 
and/or accretion surfaces. Non-thermal gyrosynchrotron radi- 
ation has been observed towards the YSO IRS 5 in the R 
CrA cloud ( Feigelson et al. 1998[ l. Synchrotron emission has 
been identified via detection of linear polarisation towards a 
very young T Tauri star by Phillips et al. (1996) (see also 
IMassi et al. 20 06 ). There is also evidence for extended syn- 
chrotron protostellar jets presented in the triple radio source in 
Serpens ( Rodriguez et al. 19"89| . 

The Corona Australis (CrA) dark cloud is a nearby (d ~ 
170 PC [ Knude & H0g"T998 f] low- and intermediate-mass star- 



forming region, where the R CrA core represent the most active 
site of star formation. The evolutionary stages range from Class 
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' There are also shorter distance estimates for the CrA star-forming 
region. Marraco & Rydgren (1981) derived a distance of only 129 pc 
for the R CrA dark cloud. Also, Casey et al. (1998) derived a similar 
distance of 129 ± 1 1 pc for TY CrA. In the present paper, the Hipparcos 
distance ~ 170 pc is adopted. 



objects to pre-main sequence (PMS ) stars ( Wilking et al. 1997 
INisini et al. 2005ilNutter et al. 2005] ). The R CrA region is there- 
fore favourable for studies of the dependence of radio continuum 
emission on the stellar evolutionary stage. 

In recent years, the cluster of YSOs in the R CrA core, also 
called the Coronet ( Taylor & Storey 1984|), have been studied 
extensively in X-rays, in the infrared, submillimetre, and in radio 
( |Hamaguchi et al. 2005t. Nutter et al. 2005.,Forbrich et al. 20061 
2007: IForbrich & Preibisch"2007l[&oppi et al. 2007| l. It has be- 
come evident that some of the radio and X-ray sources in the 
region are highly variable. Forbrich et al. (2007) performed 
the first simultaneous X-ray, radio, NIR, and optical monitor- 
ing of the cluster, and detected seven objects simultaneously in 
all wavelength bands. These series of studies have given infor- 
mation on the processes and emission mechanisms associated 
with various stages of early stellar evolution, and helped to as- 
certain the evolutionary stages of most of the cluster members. 
However, the Radio Sources 5 and 9, first detected by Brown 
(1987), have not yet been unambiguously classified. These two 
sources will be hereafter called as Brown 5 and Brown 9. 

In this paper we present centimetre continuum imaging of 
the R CrA core performed at the Australia Telescope Compact 
Array using the wavelengths 3, 6 and 20 cm. The purpose of this 
study was to determine the radio spectral indices of the various 
sources in the region, and thereby remove some uncertainty con- 
cerning their nature. At the same we add a couple of measure- 
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ment points to the time series of the varying sources. In Sect. 
2 of this paper we review briefly the radio continuum emission 
mechanisms associated with YSOs. The observations and data 
reduction procedures are described in Sect. 3. The results are 
presented in Sect. 4. In Sect. 5 we discuss the results, and in 
Sect. 6 we summarise our conclusions. 



2. Radio continuum emission from YSOs 

2.1. Spectral index and polarisation 

Radio emission from young protostars usually has a spectral 
indejfl a, in the range [-0.1, 2.0]. This range is characteris- 
tic for thermal free-free emission from partially ionized gas 
(nH^/wH -0.1, see IMartin 19961 [Anglada et al. 19 98 and ref- 
erences therein). The dependence represents optically thick 
free-free emission, whereas optically thin emission is almost 



independent of frequency (e.g..l Gudel 2002 1). Conforming with 
this emission mechanism the brightness temperatures measured 
for protostellar sources are moderate {T\, < 10"* K). The emis- 
sion can originate in ionized jets from protostars or in shocks 
produced by jets penetrating into the circumstellar material (e.g., 
ICuriel et al. 1987V Other possible sources are ionized disk wind 
(e.g., [Martin 19961 or accretion shocks on circumstellar disks 
(e.g., [Winkler & Newman"l980'; 'Ne ufeld & HoUenbach 19961 
see also [Ghavamian & Hartigan 1998) . 

The predicted spectral index of the so-called standard 
spherical stellar wind (i.e., fully ionized and isotropic stel- 
lar wind with constant velocity) is 0.6 at cm-wavelengths 
(e.g., [Panagia & Felli 1 975). According to the model of Bertout 
(1983) also accretion can produce a spectral index of 0.6. In this 
model the region emitting at radio wavelengths is photoionized 
by soft X-ray and EUV radiation from an accretion shock around 
a protostellar core. In contrast, in the model of Felli et al. (1982) 
accreting ionized circumstellar envelope produces a flat spec- 
trum with spectral indices in the range of a = -0. 1 -0.1. 

Reynolds (1986) has modelled continuum emission from 
both collimated and accelerated spherical flows. According to 
these results, a can increase above 0.6 for an accelerated flow, 
but a collimated ionized flow produces a flat spectrum with 
a < 0.6. In the model of collimated ionized flow, the opac- 
ity decreases with the distance to the central source (see also 



Anglada et al. 1998| . 



Several YSOs show negative spectral indices, a < -0.1, 
in the radio, characteristic of non-thermal gyrosynchrotron 
emission (see, e.g.. IForbrich et al. 20061 IStamatello s et al. 20071 
and references therein). Gyrosynchrotron emission arises from 
mildly relativistic electrons (Lorentz factor y < 2-3, see 
IDuIk 198 5) gyrating in magnetic fields close to the YSO and/or 
to the star-disk interaction region. Thus, non-thermal radio emis- 
sion indicates magnetic activity around protostars. 

Also shocks can accelerate electrons to mildly relativistic en- 
ergies (~ 1 MeV). Clear evidence of efficient electron accelera- 
tion associated with YSOs have been found in the Serpens triple 
source (see [Feigelson & Montmerle 1 999). Gyrosynchrotron ra- 
diation can be generated in jet induced shocks by Fermi accel- 
eration of electrons. In the so-called diffusive shock accelera- 
tion (DSA) electron scatters off" magnetic irregularities frozen 
in to the local plasma flow. This interaction can change the di- 
rection of electron, which can result in several encounters with 
the shock front enabling several encounters with the shock. 



The shock crossings lead to a systematic gain of kinetic en- 
ergy. Besides having a negative spectral index, gyrosynchrotron 
emission is characterised by rapid variability, moderate degree 
of circular polarisation {\V\/I < 20 %), and very high bright- 
ness temperature {T\, > 10''^ K) (see, e.g., [Andre et al. 19961 
IForbrich et al. 2006I I. 

Partial linear polarisation characteristic of synchrotron radi- 
ation from ultrarelativistic electrons (y » 1) was detected by 
Phillips et al. (1996) (see Sect. 1), which implies very efficient 
electron acceleration. As a distinction from the gyrosynchrotron 
radiation, the degree of circular polarisation (which to order of 
magnitude amounts to y ') of synchrotron radiation is low 
(Dulk 1985I I. 

The observed spectral index mainly describes the source av- 
erage opacity and in the case of inhomogeneous source, the ra- 
dio emission can originate from both optically thin and thick 
regions. Thus, the spectral index alone gives only a rough diag- 
nostic of the nature of YSO radio emission ([Andre et al. 19961 



lAnglada etain998i. 



2.2. Relation between radio continuum emission mechanism 
and evolutionary stage of the YSO 

The dominant radio continuum emission mechanism is sug- 
gested to change with the age of the YSO jGibb 19991 ). In the 
Class O/I stage, the principal emission mechanism is thermal 
free-free emission (accretion and outflowing jets). Continuum 
emission from Class II objects (classical T Tauri stars, CTTSs) 
at cm-wavelengths is caused by thermal free-free mechanism 
in strong ionized winds. In the Class I/II phase, possible non- 
thermal component is likely to be free-free absorbed in the sur- 
rounding dense ionized gas. The detection of a non-thermal com- 
ponent becomes more and more difficult towards higher frequen- 
cies where free-free absorption is accentuated (Giidel 2002l l. 
Later, in the Class III stage (weak-lined T Tauri stars, WTTSs), 
when circumstellar disks and ejecta are very weak or absent, 
non-thermal gyrosynchrotron emission from the exposed star 
can dominate at radio wavelengths due to increased stellar 
magnetic field (see [Feigelson & Montmerle 1999} [Gibb 19991 
Wil king et al. 200 l| l. 

Free-free emission can dominate gyrosynchrotron emis- 
sion not only if the density is high (free-free absorption), but 
also if the temperature is low or the magnetic field is weak. 
Moreover, different frequency depencies may cause switching 
between the dominant emission mechanism at a certain fre- 
quency ( Dulk 1985) . Furthermore, Gibb (1999) suggested that 
the apparent correlation between radio continuum emission and 
evolutionary stage could be due to source geometry and optical 
depth effects. 



3. Observations and data reduction 

The R CrA star-forming region was observed at 3, 6 and 20 
cm continuum with the Australia Telescope Compact Array 
(ATCAfl located near Narrabri, Australia. The observations 
were made during two observing runs in 1998 and 2000. The 
dates, configurations, phase centres, calibrators, and other ob- 
serving parameters are given in Table [1] The correlator configu- 
ration used provided at each frequency a total bandwith of 128 
MHz recorded as a 16 channel spectrum. Four products of the 
linearly polarised feeds were recorded (XX, YY, XY and YX), 



^ The spectral index is defined as 5 ^ oc v", where S is the flux density 
and V is the frequency. 



Operated by the CSIRO Australia Telescope National Facility. 



Miettinen et al.: Imaging of R CrA star-forming region 



3 



and thus all the Stokes parameters can be solved from the cali- 
brated data. 

The data were calibrated and imaged using the Miriad soft- 
ware package JSault et al. 19951 ). The visibilities were trans- 
formed into maps using natural weighting. The image decon- 
volution was done with the Miriad task mfclean. Primary beam 
correction was finally done for images (Miriad task linmos). 
Intensity maxima exceeding 3 times the local rms are considered 
as detections. For non-detections we give 3cr upper limits. 



4. Results 

The 3 and 6 cm continuum images from 1998 Jan 9/10 are shown 
in Figs. [T] and |2] respectively. The 3 and 6 cm images from 2000 
May 31 and Jun 9 toward the IRS 7 region are shown in Fig.[3j 
A comparison of the 20 and 6 cm images from 2000 Jun 10 and 
1998 Jan 9/10 showing the radio lobes associated with IRS 7B 
are presented in Fig. |4] Altogether 8 radio sources were identi- 
fied by inspecting the maps visually. The sources are listed in 
Table 121 This gives their J2000 coordinates, identifications with 
previously known sources, and suggested classifications. 

The flux densities are given in Table [3] These were obtained 
from Gaussian fits using the Miriad task imfit. For the point- 
like sources 1-6 the flux density is the peak of the Gaussian, 
whereas for the somewhat extended sources 7 and 8 we give the 
integrated intensity over the pixels above the 3cr level as calcu- 
lated by imfit. The errors are determined in the fitting proce- 
dure. As can be seen from Tabled only Brown 5, IRS 5, and the 
IRS 7A/B pair are visible at all three wavelengths. 

We also obtained spectral indices between 3 and 6 cm 
(a?'^™), and 6 and 20 cm (q'3°'^™). For the determination of these 

^ 3cm^' ^ 6cm ' 

spectral indices, the maps at the two frequencies were convolved 
to the same angular resolution using the Miriad task convol. 
On May 3 1 and Jun 9, 2000 we used configurations which gave 
similar synthesized beams at 3 and 6 cm in order to derive re- 
liable spectral indices for extended sources. In 1998, Jan 9/10, 
the 3 and 6 cm observations were done simultaneously, and in 
this case the spectral indices should be more reliable for variable 
point sources. The spectral indices Q'3™ and ag*^™, and their Icr 
errors are given in Table |4] The error estimates are based on the 
rms errors of the flux densities, and do not take into account the 
possible errors caused by different beams or source variability. 



5. Discussion 

All the eight sources detected in the present study are known 
from previous radio surveys ( IBrown 19871 (2 and 6 cm, VLA); 
ISuters et al. 19961 (2, 3.5, 6, and 20 cm, VLA and ATCA); 
[Feigelso neTal. 1998| (3.5 cm, VLA) ; IChoi & Tatematsu 20041 
(0.7 cm, VLA); IForbrich et al.^006t 2007 (3.6 cm, VLA)). 
Five of the sources, IRS 5, IRS 6, IRS 1, R CrA, and IRS 
7A, have counterparts both in the near/mid-infrared and X- 
rays. Moreover, IRS 7B, has been recently detected in the mid- 
infrared (but not in the NIR) and X-rays. Brown 9 corresponds 
to a submm peak ( [Groppi et al. 2007) . Brown 5 has only been 
detected in radio, but it may be associated with a near-IR source 
( WiUdng et al. 1997 1. In what follows we discuss the radio prop- 
erties of the detected sources in detail. 
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Fig. 1. Map of the 3 cm continuum made from the 1998 Jan 9/10 
visibility data toward the R CrA star-forming region. The loga- 
rithmic grey-scale range from 0.15 to 5 mJy beam"'. The con- 
tours are from 0.29 mJy beam"' (3cr) to 2.9 mJy beam"' in steps 
of 0.29 mJy beam"'. Shown at the bottom right-hand corner is 
the synthesized beam: FWHM = 3'.'4 x 3'.'1 and RA.= -31°.0. 
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Fig. 2. Map of the 6 cm continuum made from the 1998 Jan 9/10 
visibility data toward the R CrA star-forming region. The loga- 
rithmic grey-scale range from 0. 1 to 5 mJy beam"' . The contours 
are from 0.24 mJy beam"' (3cr) to 2.4 mJy beam"' in steps of 
0.24 mJy beam"'. Shown at the bottom right-hand corner is the 
synthesized beam: FWHM = 6'.'5 x 5'.'7 and RA.= -19°6. 



5. 1 . Comments on individual sources 
5.1.1. The IRS 7 region 

IRS 7 is the most active complex in the Coronet cluster A 
dense, bipolar molecular outflow and a large rotating molecu- 
lar disk (outer radius > 3900 AU at li = 170 pc) were detected 
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Date of observations 1998 Jan 9/10 2000 May 28 2000 May 31 2000 Jun 9 2000JunlO 

Time (UT) intei-val 19:06:05-08:40:55 09:56:35-23:11:45 09:33:35-21:56:25 09:26:15-22:21:55 09:25:55-22:03:35 

Wavelength 3 and 6 cm 3 cm 3 cm 6 cm 20 cm 

Bandwith 128 (16x8) MHz;^ 

Configuration 1.5A ~ 1.5A 1.5A 6B 6B 

Hi'-distance range 5-149 and 2.6-74.5 U 5-149 U 5-149 U 3.6-99.5 U 1-30 M 

Phase centre i9hoi">in44M 19''01™»56?5 19''0r™56^5 19''0r™5ffi5 19''01"™56?5 

(0-2000.0,^2000.0) -36°58'40'.'3 -36°57'28" -36°57'28" -36°57'28" -36°57'28" 

Flux calibratoi-* PKS 1934-638 PKS 1934-638 PKS 1934-638 PKS 1934-638 PKS 1934-638 

Phase calibrator? B1849-36 B1933-400 B1933-400 B1933-400 B1933-400 

Image processing Naturally weighted" 

FWHM of synthesized beam 3'.' 4 X 3'.' 1 and 6'.' 5 X 5'.' 7 4'.'4 X 3'.'3 4'.'4 X 3'.'2 4'.'2 X 2'.'3 13'.'lx6'.'6 

Beam position angle -31?0 and -19^6^ 2?5 1?9 -2?8 0;i 

rms noise level [//Jy beam"'] 97and81 76 47 44 120 



" Used in all cases. 

* Calibrator flux densities are the following: PKS 1934-638: Sscm = 2.84 Jy, Sgcm = 5.83 Jy, ^zocm = 14.95 Jy; B1849-36: 
Sscm = 0.42 Jy, 56cm = 0.70 Jy; B1933-400: 53cm = 2.09 Jy, 56cm = 1-38 Jy, 5 20cm = 0.81 Jy (from the ATNF web site 
|htt p : //www .narrabri . atnf . csiro . au/calibrators/ I. 

The difference between the beam position angles is due to the different amount of flagged uv data. 



Table 2. Positions, identification and the class of the detected radio sources. 



Peak position^ 



No. 


«20oo.o [h:m:s] 


^2000.0 [":':"] 


Source ID 


Source class 


1 


19 01 43.3 


-36 59 13 


Brown 5_^ 


unclear 


2 


19 01 48.0 


-36 57 22 


IRS5~ 


Class I 


3 


19 01 50.5 


-36 56 38 


IRS 6 


Class II (CTTS) 


4 


19 01 50.7 


-36 58 10 


IRS 1 


Class I 


5 


19 01 53.6 


-36 57 12 


RCrA 


HAe 


6 


19 01 55.2 


-36 57 17 


Brown 9* 


Class ? 


7 


19 01 55.3 


-36 57 22 


IRS 7A, IRS 7W 


Class I 


8 


19 01 56.4 


-36 57 28 


IRS 7B, IRS 7E~ 


Class O/I 



" These positions refers to peak positions at 3 cm observed in January 1998. 
* Source number in Brown 1987 
Source used in .Forbrich et al. 2006i 2007. 

Table 3. Flux densities with +lcr eiTors of the detected radio sources at 3, 6 and 20 cm. The values are corrected for the primary 
beam response. 



Source 


1998 Jan 9/10 


2000 May 28 


2000 May 31 


2000 Jun 9 


2000 Jun 10 




53 cm [mjy] 


56 cm [mJy] 


53 cm [mJy] 


53 cm [mJy] 


56 cm [mJy] 


520 cm [mJy] 


Brown 5 


0.96 ± 0.06 


1.42 ±0.10 


1.40 ±0.13 


1.49 ±0.20 


1.05 ±0.01 


1.80 ±0.02 


IRS 5 


0.81 ±0.03 


1.06 ±0.04 


1.87 ±0.21 


0.70 ± 0.04 


0.62 ± 0.03 


0.40 ± 0.04 


IRS 6 




0.46 ± 0.04 










IRS 1 


0.33 + 0.04 




0.62 ±0.21 


0.52 ± 0.06 


0.21 ±0.02 




RCrA 


0.39 ± 0.03 




0.46 ± 0.02 


0.35 ± 0.04 


0.29 ± 0.03 




Brown 9_^ 


0.74 ±0.14 




1.88 ±0.08 


1.93 ±0.10 


1.58±0.11 




IRS 7a5 


4.19 ±0.07 


5.49 ±0.12 


6.00 ± 0.02 


5.90 ±0.01 


5.60 ± 0.09 


5.10 ±0.08 


IRS 7B~ 


2.56 ± 0.06 


2.75 ± 0.03 


3.65 ± 0.07 


3.70 ± 0.07 


3.42 ± 0.06 


3.67 ± 0.06 



° Blended with IRS 7A at 6 cm in the dataset 1998 Jan 9/10, and at 20 cm in the dataset 2000 Jun 10. 

* For these extended sources we give the integrated flux densities over the pixels above the 3cr level as calculated by the Miriad task imfit. 



in the IRS 7 region (I Anderson et al. 1997al 1997b). Groppi et 
al. (2004) found that the envelope surrounding the disk show 
signs of infall in the HCO^(4 - 3) spectra. The region contains 
three compact radio sources, IRS 7A, IRS 7B, and Brown 9. 
The rotation axis of the molecular envelope coincides with a 
string of HH objects (Hartigan & Graham 1987; see Fig. 10 in 
[Anderson et al. 19 97a). These features suggest ongoing power- 
ful accretion and possible shocks associated with infall onto an 
accretion disk. 

IRS 7A and IRS 7B are the brightest radio sources in the 
whole cloud. The evolutionary stages of these sources have been 



discussed by Wilking et al. (1997), Choi & Tatematsu (2004), 
Hamaguchi et al. (2005), Forbrich et al. (2006), and Groppi et 
al. (2007). Brown 9 (IBrown 1987 ) is located 5" north of IRS 7A 
(850 AU at = 170 pc). Brown 9 and IRS 7A are blended in the 
6 cm map from 1998 and in the 20 cm map from 2000. 

IRS 7A has not shown much radio variability. The flux den- 
sities of IRS 7A are very similar to those observed by Brown 
(1987), Feigelson et al. (1998) and Forbrich et al. (2006). The 
flat spectrum with = 0.19 + 0.04 and 0.11 + 0.03 and 

ffg^™ = 0.34 + 0.02 can be explained by optically thin emis- 
sion from ionized circumstellar envelope or a collimated jet. IRS 
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Table 4. Spectral indices with +lcr errors of the observed sources. 



Source 


^6 cm (Y 


^6 cm b 
^cm™ 


^20 cm c 


Brown 5 


-0.34 + 0.12 


0.76 ± 0.29 


-0.30 d 


:0.01 


IRS 5 


-0.12 ±0.07 


0.26 ±0.16 


0.54 ± 


0.08 


IRS 1 




1.96 ±0.32 






RCrA 




0.41 ±0.33 






Brown 9 




0.43 ±0.19 






IRS 7A 


0.19 ±0.04 


0.11 ±0.03 


0.34 ± 


0.02 


IRS 7B 


0.38 ± 0.02 


0.17 ±0.06 


-0.05 d 


:0.01 



" ffj™ derived from 1998 Jan 9/10 data using convolved flux densities. 
* '^icm derived from 2000 May 31 & Jun 9 data using original flux densities, 
ttj^™ derived from 2000 Jun 9/10 data using convolved flux densities. 



7A has been detected both in the near-infrared and in the X- 
ray band. In view of the relatively bright radio emission, these 
characteristics suggest that the source represents the protostellar 
Class I. 

The relationship between IRS 7A, Brown 9 and the nearby 
10//m source between the two ( Wilking et al. 1997 1 has been dis- 



cussed in Choi & Tatematsu (2004, hereafter CT). They suggest 
that these sources belong to a system formed by a central star 
(close to the lOyum source, associated with the A = 6.9 mm 
source No. 3 of CT) and ionized jets (IRS 7A and Brown 9) on 
its both sides. By inspecting the map of CT, the A = 3.6 cm map 
of Forbrich et al. (2006), and the SCUBA 450 jum continuum 
map of Groppi et al. (2007), one agrees that the three sources 
can be physically connected. It is, however, unlikely that the X- 
ray source IRS 7A and the submm source Brown 9 would just 
be jets arising from the lOyum source associated to IRS 7. In all 
other near-IR sources of our sample the hottest spot, i.e. the X- 
ray source, shows the location of the central protostar. On the 
other hand, the submm/far-IR spectrum of Brown 9 is not possi- 
ble to explain with free-free emission from plasma, but it is char- 
acteristic of a Class protostar (see Fig. 8 in Groppi et al. 2007 1. 

The flux densities of Brown 9 at 3 cm found in the present 
study are mostly similar to those found in the previous stud- 
ies ( [Feigelson et al. 1998} IForbrich et al. 20061) . However, the 
source shows a significant brightening (by a factor of two) since 
the observation of Brown (1987). The spectral index, Q'3™ - 
0.43 + 0.19, is consistent with thermal free-free emission orig- 
inating in a collimated ionized flow. The source of this flow is 
likely to be embedded in core detected with SCUBA. 

The IRS 7B flux densities derived here are ~0.5-1.5 mJy 
higher than those reported by Brown (1987), Feigelson et al. 
(1998) and Forbrich et al. (2006). This difference can be at- 
tributed to source variablity. The spectral indices between 3 and 
6 cm are a^™ = 0.38 + 0.02 and 0.17 + 0.06. The spectral index 
between 6 and 20 cm, a?"™ = -0.05 + 0.01, is close to the limit 

' o cm ' 

of thermal emission and non-thermal gyrosynchrotron emission 
(a - -0.1). We note, however, that c?'^'" is derived from simul- 

^ ^ ' ' 3 cm 

taneous observations (1998, Jan 9/10) and is much more signif- 
icant than the value for ag^™ (5cr) based on observations from 
two consecutive days. 

IRS 7B is associated with weak, extended radio emission. 
The lobes on the northeastern and southern sides of the point 
source are clearly visible on the 6 cm map from 1998 and the 20 
cm map from 2000 (see Fig.|4|l. These two maps were obtained 
using configurations with the shortest Mv-distances (see Table[T]). 
Traces of these lobes can be seen in the 3 and 6 cm images from 
2000 (Fig. O, which represent the highest angular resolution 
available here. The lobes become visible also in these data by 



applying a Gaussian taper which produces a resolution similar 
to that in the 20 cm image. The S/N ratio is, however, too poor 
to compare individual features at 3 and 6 cm. The extended emis- 
sion was already detected by Brown (1987) at 6 cm and by Suters 
et al. (1996) at 20 cm. The radio lobes may be associated with a 
stage prior to the formation of a HH object (Curiel e t al. 1993T l. 
The fact that the lobes are more prominent at 6 and 20 cm than at 
3 cm or shorter wavelengths suggests a negative spectral index 
and consequenly, a possible gyrosynchrotron emission compo- 
nent. This led Harju et al. (2001) to suggest that the source could 
be a distant radio galaxy or a Galactic microquasar. Later ob- 
servations have, however, ruled out this possibility. Pontoppidan 
et al. (2003) detected IRS 7B in the 4.67 jjm absorption lines 
of solid CO, which clearly originate in circumstellar material. 
Moreover, Hamaguchi et al. (2005) pointed out that it would be 
very unlike to observe AGN with X-ray flux of IRS 7B in the 
IRS 7 region. Also the X-ray spectra of IRS 7B does not resem- 
ble those of AGNs. Many properties of IRS 7B, e.g. the large 
equivalent width of the fluorescent X-ray iron line and the spec- 
tral energy distribution, suggest it to be a Class source or an ob- 



ject between Class and Class I phase (Hamaguchi et al. 2005 
INutter et al. 20051 [Groppi et al. 2007| l. 



Given the spectral indices and other properties of IRS 7B, 
the continuum radiation is probably coming from the surround- 
ing gas which is shock-ionized by the high velocity jet. The ra- 
dio lobes associated with the thermal radio jet possible have a 
gyrosynchrotron emission component. This could be understood 
in terms of DSA of electrons (see Sect. 2.1). 



5.1 .2. Brown 5 - a brown dwarf ? 

This is the Radio Source 5 detected by Brown (1987) using the 
VLA at 6 cm. It has been suggested that the infrared coun- 
terpart of Brown 5 could be a brown dwarf if it lies in the R 
CrA cloud (Wilking et al. 1997HFeigelson et al. l"998t see also 
Forbrich et al. 20061. The infrared source has been detected in 
the H and K bands only, with K = 16.4, and H - K = 1.45. The 
value of H - K suggests a large interstellar or circumstellar red- 
dening ((H - K}o ~ 0.5 for an M9V dwarf, [Wilking et al. 1999| l. 
Therefore the object is likely to lie in the cloud or behind it. It 
should be noted that brown dwarfs are seldom detected in radio 
(see, e.g., Giidel 2002. [Osten & Jayawardhana 2006) . During 
flares the flux densities of nearby (d ~ 10 pc) brown dwarfs can 
increase to a few millijanskys, but the average value is an order 
of magnitude lower (radio luminosity Lv,5ghz < 10'^ ergs s"' 
Hz"'). Assuming that the object lies at the distance of the CrA 
cloud its (nearly constant) radio luminosity is Lv.sghz ~ 3 ■ 10'^ 
ergs s"' Hz"'. This high luminosity would be characteristic of a 
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Fig. 3. Top: The 3 cm contour map made from the 2000 May 3 1 
visibility data toward the IRS 7 region. The contours are from 
0.14 mJy beam"' (3cr) to 1.4 mJy beam"' in steps of 0.14 mJy 
beam"'. The logarithmic grey-scale range from 0.15 to 5 mJy 
beam"'. Shown at the bottom right-hand corner is the synthe- 
sized beam: FWHM = 4'.'4 x 3'.'2 and P.A.= 1°9. Bottom: The 6 
cm contour map made from the 2000 Jun 9 visibility data toward 
the IRS 7 region. The contours are from 0.13 mJy beam"' (3cr) 
to 1.3 mJy beam"' in steps of0.13mJy beam"'. The logarithmic 
grey-scale range from 0.1 to 5 mJy beam"' . Shown at the bottom 
right-hand corner is the synthesized beam: FWHM = 4'.'2 x 2'.'3 
and RA.= -2?8. 



magnetically active star rather than a brown dwarf ("Giide l 20021) . 
On the other hand, the fact that Brown 5 has not been detected in 
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Fig. 4. Top: The 20 cm contour map made from the 2000 Jun 
10 visibility data toward the IRS 7 region. The contours are 
from 0.36 mJy beam"' (3cr) to 3.6 mJy beam"' in steps of 0.36 
5mJy beam"'. The logarithmic grey-scale range from 0.23 to 7 
mJy beam"'. White crosses shows the 3 and 6 cm peak posi- 
tions. Shown at the bottom right-hand corner is the synthesized 
beam: FWHM = 13'.'1 x 6'.'6 and RA.= 0°1. Bottom: Blow-up 
of Fig. |2] showing the IRS 7 region. Shown at the bottom right- 
hand corner is the synthesized beam: FWHM - 6'.'5 x 5'.'7 and 
RA.= -19%. 



X-rays dForbrich et al. 20061 IForbrich et al. 2007] l makes coro- 
nal emission from a late-type star a less likely alternative. 

The fluxes observed here are similar at all wavelengths to 
those observed by Suters et al. (1996), Feigelson et al. (1998) 
and Forbrich et al. (2006). On one day Suters et al. detected an 



Miettinen et al.: Imaging of R CrA star-forming region 



7 



unusually high 3 cm flux, and later on, Forbrich et al. (2006) 
found that the source is variable on timescales of days to months. 



The radio spectral indices, a"""" = -0.34 + 0.12 and orf ™ = 

r '3 cm 6 cm 

-0.30 + 0.01, indicate optically thin gyrosynchrotron emission, 
and so magnetic field activity in this object. The former spectral 
index is similar to that found by Suters et al. on Jul 12, 1992 
(-0.28 + 0.30), but our 6/20 cm index is very different from their 
values (af™ -0.1). Spectral index derived from 2000 May 3 1 / 
Jun 9 data, a^™ = 0.76 + 0.29, is probably not meaningful since 
the source is variable with a scale of days. 

According to Suters et al. (1996) the high obscuration 
and rapid variations make it more likely that Brown 5 is a 
PMS star rather than an extragalactic source. However, extra- 
galactic nature cannot be totally ruled out as interstellar scin- 
tillation can cause so called IntraDay Variability (IDV, see 



Gabanyietal. 2007 1. 



5.1.3. IRS 5 

The Class I protostar IRS 5 is known to be a highly variable radio 
source. Our ATCA observations yield mainly similar flux den- 
sities to those observed by Brown (1987), Suters et al. (1996), 
Feigelson et al. (1998) and Forbrich et al. (2006; 2007). At 20 
cm we found a flux of 0.40 mJy, over 2 times lower than the 
lowest value found by Suters et al. (1996). Suters et al. detected 
a flare of 6.6 mJy at 6 cm in 1992 and concluded that IRS 5 
experienced an outburts in mass loss. The radio spectral indices 
of IRS 5 found in the present study are = -0.12 + 0.07 

and 0.26 + 0.16 and af"'" = 0.54 + 0.08. Wenote that the latter 

6 cm 

two values can be deceptive since the source is highly variable. 
Suters et al. found spectral indices Q'3™ ~ 0.1 and ~ 0.9, and 
20 cm varied from ~ -0.01 to -0.1, all very different from those 

6 cm ' 

derived here. On the other hand, the indices of Suters et al. have 
rather large errors. 

IRS 5 has been detected in X-rays (Koyama et al. 
1996; INeuhauser & Preibisch 19971 ), and recently Forbrich et 
al. (2006; 2007) found it to be the most variable source in X- 
rays in the R CrA region (see also Forbrich & Preibisch 20071 ). 
Feigelson et al. (1998) detected circularly polarised non-thermal 
emission from IRS 5 at 3.5 cm, which is also the first such de- 
tection for a YSO. This was later confirmed by Forbrich et al. 
(2006), but interestingly, there was also occasion when IRS 5 
did not show circular polarisation in their observations. In the 
present study, the Stokes-V emission of IRS 5 was detected only 
at 3 cm on May 31, 2000, with the flux density of 0.23 + 0.04 
mJy, the polarisation fraction, V/I, being as large as about 33 %. 
The 3cr upper limits of V emission at 3 and 6 cm observed in 
1998 Jan 9/10, at 3 cm in 2000 May 28, and at 6 and 20 cm in 
2000 Jun 9/10 are about 0.07, 0.20, 0.18, 0.06, and 0.04 mJy, 
respectively. The corresponding upper limits of V/I are about 9 
and 19 % at 3 and 6 cm in 1998 Jan 9/10, respectively, and about 
10 % during other measurements in 2000. 

The slightly negative spectral index between 3 and 6 cm, 
and circular polarisation suggest that the radio emission has a 
substantial non-thermal gyrosynchrotron component. This con- 
forms with the fact that IRS 5 does not appear to be power- 
ing an outflow (Feigels on et al. 1998[ see also |Wang et al. 2004] l. 
Detectable non-thermal emission implies that the immediate 
surroundings of the star are relatively free of optically thick 
ionized gas (see Sect. 2.2). The classification of IRS 5 is 
complicated by the fact that it is a binary with a separa- 
tion of ~ (y.'6 or 102 AU at 170 pc dChen & Gra ham 1993; 
INisini et al. 20051 ). The binary is also marginally resolved in X- 



rays jForbrich & Preibisch 20071 ). This together with the rapid 
variability make the determination of the spectral indices uncer- 
tain. Nevertheless, the radio properties of this source suggest that 
at least one of the components has advanced beyond the Class I 
stage. 



5.1 .4. Other detections 

The T Tauri star IRS 6 was detected only at 6 cm in the 1998 
visibility data with the flux density 0.46 + 0.04 mJy. Thus, the 
spectral index cannot be determined. The 3cr upper limits of the 
3 cm flux density ai-e 0.02 and 0.03 mJy (on Jan 9/10, 1998, and 
May 28/31, 2000, respectively), and the 3cr upper limits of the 
6 and 20 cm flux densities are 0.13 and 0.10 mJy (on Jun 9/10, 
2000). There are no reported radio emission from this source 
before Forbrich et al. (2006). Forbrich et al. (2006; 2007) also 
found very low fluxes of ~ 0.1 mJy at 3.6 cm. 

IRS 1 is another Class I object in the Coronet cluster 
( |Wilking et al. 1986[[T992b . It is also known as HH 100-IR, and 
is Ukely a driving source of the HH flow consisting of HH 96, 
97, 100, and 98. The flux densities of IRS 1 in the present 
work are very similar to those observed by Brown (1987), 
Suters et al. (1996), Feigelson et al. (1998) and Forbrich et 
al. (2006). At 20 cm IRS 1 remained below the 3o- level of 



0.36 mJy. The spectral index ffj™ = 1.96 ± 0.32 is consis- 
tent with the value - 1.67 + 0.72 derived by Suters et al. 

3 cm 

(1996). This high spectral index suggests optically thick free- 
free emission, which can attributed to shocked material. IRS 1 
has been detected as a hard X-ray source (Koyama et al. 1996; 
Neuhauser & Preibisch 199"7l ). Even though IRS 1 does not show 
significant radio variabiUty, it shows clear signs of X-ray vari- 
ability (I Suters et al. 19961 IForbrich et al. 2006i According to 
Nisini et al. (2005) a very high fraction of the bolometric lu- 
minosity of IRS 1 is due to accretion. In view of these results it 
seems possible that the cm emission originates in an accretion 
shock in the circumstellar disk. 

The Herbig Ae star R CrA has remained undetected 
in several surveys, indicating long-term variability (see 
IForb ric h et al. 20061 and references therein). Our 3 cm flux den- 
sities are somewhat higher than those observed by Forbrich et 
al. (2006; 2007). The 6 cm flux of R CrA was 0.29 + 0.03 mJy. 
The source was not detected at 20 cm above the 3cr level of 0.36 
mJy. 

Koyama et al. (1996) detected hard X-ray emission close 
to the position of R CrA and Forbrich et al. (2006) clearly de- 
tected X-ray emission from R CrA in all of their Chandra and 
XMM-Newton datasets. As Forbrich et al. (2006) pointed out, 
the X-ray emission from this star does not fit into our under- 
standing of stellar X-ray emission mechanisms (magnetically 
driven corona or shocks caused by strong stellar wind; see also 
IForbrich & Preibisch 20"07l ). 

Given the substantial error in the spectral index, Q'3™ = 
0.41 + 0.33, it is not possible to make firm conclusions about 
the emission mechanism for this continuum radiation. However, 
since R CrA should not drive very strong stellar wind, the up- 
per limit of the spectral index (~ 0.7) is excluded. On the other 
hand, since R CrA is likely to be a driving source of coUimated 
outflow (see, e.g., ISuters et al. 19961 ). radio emission is probably 
produced by a coUimated, ionized flow. 
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5.1.5. Non-detections 

IRS 2: The Class I protostar IRS 2 is one of the strongest X-ray 
sources among the objects studied by Forbrich et al. (2006). The 
non-detection is probably due to primary beam attenuation far 
from the observation's phase centre. 

IRS 9: The Class I protostar IRS 9 remains undetected at ra- 
dio also in previous radio observations (see lForbrich et al. 20061) . 

T CrA: The Herbig Ae star, T CrA, of spectral FOe is one of 
the most massive stars in the R CrA core. The non-detection is 
consistent with the earlier studies of Brown (1987) and Suters et 
al. (1996). 



X-rays to radio, are required to understand the nature of the 
youngest protostars (IForbrich et al. 20071 Groppi et al. 2007 1. 
At the same it is evident that in the radio regime, higher sen- 
sitivity and angular resolution than available now are needed to 
detect a significant fraction of YSOs in a dense cluster like R 
CrA, and to study the physical processes which determine their 
evolution. 
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6. Summary and conclusions 

We have conducted a survey of radio continuum emission at 3, 
6 and 20 cm in the R CrA star-forming region using the ATCA. 
Altogether eight pointlike sources were detected. Seven of these 
can be assigned to YSOs, having counterparts in the X-rays, in- 
frared or submm. We studied their spectral indices, variability, 
and polarisation to gain information on possible emission mech- 
anisms. In this manner we examined the classifications of the 
YSOs in the light of the evolutionary scenario presented by Gibb 
et al. (1999), which relates radio properties to protostellar stages. 

Probably the youngest protostars in our sample are Brown 9 
and IRS 7B. Brown 9 is a Class candidate and IRS 7B has been 
suggested to represent a transition stage between Class and 
Class I ( Groppi etal. 2007 1. The spectral indices of these two 
sources are consistent with free-free emission from a coUimated 
jet. IRS 7B shows signs of extended lobes at 6 and 20 cm. The 
lobes may have a gyrosynchrotron component which possibly 
can be understood in terms of Fermi acceleration of electrons in 
shocks associated with high velocity jets from the protostar 

The three Class I protostars detected in this survey, IRS 7A, 
IRS 5 and IRS 1 have very different characteristics. The spec- 
tra of IRS 1 and IRS 7A are consistent with free-free emission, 
with the distinction that emission from IRS 1 is optically thick, 
possibly originating in shocks, whereas emission from IRS 7A is 
optically thin, suggesting a collimated jet. Radio emission from 
the wide binary star IRS 5 has a substantial non-thermal compo- 
nent indicated by a negative spectral index, rapid variability, and 
strong Stokes-V emission on one of days of observation. These 
characteristics are signs of magnetic activity close to the stellar 
photosphere. This suggests that at least one of the components 
of IRS 5 is free of absorbing ionized gas, and may have akeady 
reached an advanced stage. 

Two infared-bright Class I candidates in the region, IRS 2 
and IRS 9, remained undetected in radio. In accordance with pre- 
vious results, R CrA (HAe) and the Class II star IRS 6 (CTTS) 
show weak radio emission. The so far unclassified bright radio 
source Brown 5 has a mildly negative spectral index, character- 
istic of optically thin gyrosynchrotron radiation. This source is 
likely to be an obscured Class III star or an extragalactic source. 

The present observations are in rough agreement with the 
scheme of Gibb et al. (1999) where the dominant emission 
mechanism changes with the age of an YSO. It was found, how- 
ever, that Class I stars form a very heterogeneous group as to 
their radio properties. This can be perhaps understood in terms 
of a rapid decline in the jet activity, which also changes the effi- 
ciency of free-free absorption. The results from IRS 7B suggest 
that energetic jets from very young protostars can give rise to 
non-thermal emission via shock acceleration. 

The experiences from the R CrA stellar cluster have shown 
that observations in a wide range of wavelengths, in fact from 
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